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REAL-TIME IMAGING SYSTEM AND METHOD 

Applicants claim the benefit of U.S. Provisional Application No. 60/310,082 filed 
August 7, 2001, the entire disclosure of which is hereby incorporated herein. 

Field of the Invention 

The present invention is directed to a real-time imaging system and method that is 
particularly useful in the medical field, and more particularly, to a system and method for 
imaging and analysis of tissue using optical coherence tomography. 

Bacl^round of the Invention 

A variety of imaging techniques are used for the medical diagnosis and treatnient of 
patients. Ultrasound imaging represents a prevalent teclmique. Ultrasound uses soxmd waves to 
obtain a cross-sectional image of an object. These waves are radiated by a transducer, directed 
into the tissues of a patient, and reflected fi-om the tissues. The transducer also operates as a 
receiver to receive the reflected waves and electronically process them for ultimate display. 

Another imaging teclmique is referred to as Optical Coherence Tomography (OCT). OCT 
uses light to obtain a cross-sectional image of tissue. The use of light allows for faster scanning 
times than occurs in ultrasoimd technology. The depth of tissue scan in OCT is based on low 
coherence interferometry. Low coherence inteiferometry involves splitting a hght beam firom a 
low coherence Ught source into two beams, a sampling beam and a reference beam. These two 
beams are then used to form an interferometer. The sampling beam hits and penetrates the tissue, 
or other object, under measurement. The sampling or measurement beam is reflected or scattered 
firom the tissue, carrying information about the reflecting points fix)m the surface and the depth of 
fesu e. The r ef^nGe4?eam-hite-^-ef^^3ce^Fefl ector » s uch a s» f o r e x am pter-a- mirror - or a 
dififraction grating, and reflects from the refiarence reflector. The reference reflector either moves 
or is designed such that the reflection occurs at different distances from the beam splitting point 
and returns at a different point in time or in space, which actually represents the depth of scan. 
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The time for the reference beam to return represents the desirable depth of penetration of tissue 
by the sampling beam. 

When the reflected beams meet, mtensities j&om respective points with equal tune delay 
foiTO interference. A photodetector detects this interference and converts it into electrical signals. 
The signals are electronically processed and ultimately displayed, for example, on a computer 
screCTi or other monitor. 

Summary of the invention 

The present invention provides a system and method for overcoming or minimizing the 
problems of prior optical coherence tomography systems and improving on other imaging 
methodologies. Software techniques that are used for real-time imaging in OCT (Optical 
coherence tomography), particularly for correcting geometric and angular image distortions, hi 
addition, a methodology for quantitative image correction in OCT images mcludes procedures 
for correction of non-telocentric scan patterns, as well as a novel ^roach for refraction 
correction in layered media based on Fennat's principle. 

The foregoing and other features of the invention are hereinafter fully described and 
particularly pointed out in the claims, the foUowmg description and annexed drawings setting 
forth in detail certain illustrative embodiments of the invention, these embodiments being indica- 
tive, however, of but a few of the various ways in which the prmciples of the invention may be 
employed. 

Brief Description of the Drawings 
FIG. 1 is a timing diagram for double-sided line acquisition. 
•^qEG. 2 i s a mappiag-amy-used^n-badgwafd-tet^ 

FIG. 3 is an example of inserted zoom realized with mapping arrays in real-time. 
FIG. 4 is a flow chart for tlie determination of whether a pixel is displayed as a structural 
or flow value. 

FIG. 5 illustrates an exemplary non-linear scanning. 
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FIG. 6. 1 is an OCT system in accordance with the present invention illustrating 
components and synchronization. The thick Unes represent optical signals, dash Unes represent 
electronic signals, and thin Unes represent synchronization signals. 

FIG. 6.2 schematically illustrates optical power-conserving interferometer configurations. 

FIG. 6.3 is a timing diagram for OCT synchronization electronics. 

FIG. 6.4 is a block diagram of an endoscopic OCT (EOCT) system. Light from a high- 
power broadband source is coupled through an optical circulator to a fiber-optic Michelson 
interferometer. The EOCT catheter probe and probe control imit constitute one arm of the 
inlerferometer, and a rapid-scanning optical delay Une constitutes the other arm. The gray lines 
represent optical paths and black lines represent electronic paths. 

FIG. 6.5 is a block diagram of a simple frame grabber. Video signals can be either 
composite or non-composite. External sync signals are selected by tlie acquisition/window 
control circuitry. 

FIG. 6.6 signal components (normal) and input controls (italics) of the horizontal video 

signal. 

FIG. 6.7 signal components (normal) and input controls (italics) of the vertical video 

signal. 

FIG. 6.8 plot of the retinal false-color scale represented in RGB color space. Green and 
blue color values are identical between 209-255 pixel values. 

FIG. 6.9 comparison of an in vivo human retinal OCT image along the papillomacular 
axis represented in a) linear grayscale, b) inverse linear grayscale, and c) retinal false-color scale 
(with labels). 

FIG. 6.10 is a one-dimensional forward and inverse mappings. 
^4G7-67l44^11ustrat^s4mage-rert;ation-^ransfbm 
FIG. 6.12 illustrates rectangular to polar coordinate transformation. 
FIG. 6.13 is a timing diagram for double-sided line acquisition. 
FIG. 6.14 is a motion artifact reduction by cross-correlation scan registration. Patient 
axial motion during the original retinal image (a) acquisition was estimated from a profile built 
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up from the peak of the cross-correlations of each A-scan with respect its neighbor (b). The 
resulting profile is then hi^-pass filtered to preserve real retinal profile, and used to re-register 
each individual A-scan in the image (c). 

FIG. 6.15 is a schematic of systems theory model for OCT. 

FIG. 6.16 is an example of digital deconvolution of low-coherence interferometry data. 
FIG. 6.16a is an observed interferogram of a cover slip resting on a piece of paper. 
FIG. 6. 16b is an interferogram obtained with a mirror in the sample arm. 
FIG. 6.16g is a deconvolved impulse response profile. 

FIG. 6.17 illustrates original (a) magnitude-only deconvolved (b) and iteratively 
deconvoled (c) cross-sectional OCT images revealing cellular structure in a onion sample. Both 
deconvolutions resulting in a resolution increased by a factor of approximately 1 .5 or 
approximately 8 micrometers FWHM resolution in the deconvolved images, although the 
iterative restoration algorithm preserved image dynamic range significantly better. 

FIG. 6.18 is an illustration of coherent OCT deconvolution. In the top image, the 
magnitude and phase of a demodulated OCT A-scan data obtained from two closely spaced 
glass-air interfaces with slightly distinct separation, results in (a) destractive (note 180'' phase 
shift at the mid-point) and (b) constructive interference between the reflections. In the middle 
images, deconvolution of the A-scan data performed using only the magnitude of the input data, 
leads to inaccurate positioning of reflections and spurious reflections in flie calculated impulse 
response. In the bottom image, complex deconvolution recovers the true locations of the 
intCTfaces in both cases and thus enhances resolution by a factor of approximately 1 .5, as well as 
reducing speckle noise. 

FIG. 6.19 is a demonstration of depth-resolved OCT spectroscopy in a discrete optical 
element. The spectral transfer characteristic |H(k)|^ of the Ught reflected from 1) the front 
-surfaee'and^)-the-rear-siHface-of-^-^^ 

were obtained by digital processing of windowed OCT A-scans of the filter. The measured 
spectral widths correspond well with the manufecturer's specifications (SLD spectral width 
47 nm FWHM; filter bandwidth nm FWHM). 

FIG. 6.20 is a table of usefiil spatial transformation (point-set operation) matrices. 
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FIG. 7 is an illustration of using a pointer array as a mapping array to allow for fast 
backward transformation. 

FIG. 8 is an illustration for correction for sinusoidal scanning. 
FIG. 9 illustrates correction for divergence. 

FIG. 9a indicates coordinates and measures in the intermediate image, and b) provides the 
same for the target image. 

FIG. 10 illustrates a path of light through different layers of tissue, refracted at the points 
Pbi and Pb2. L = Li H- L2 + L3 should be minimized to JBnd the right path. 

FIG. 1 1 includes a series of OCT images of the temporal anterior chamber angle of a 
human eye, imaged in vivo at 8 §)s, in different stages of dewarping. 

FIG- 1 la is a raw image. 

FIG. 1 lb illustrates removal of nonlinear reference mirror movement, FIG. 1 Ic illustrates 
divergence correction of a handheld scanner, FIG. lid correction for refraction at the air-comea 
boundary (n^omca = 1-38), FIG. lie additionally corrects for refraction at the endotheUum-aqueous 
boundary (nacqueous - 13 3). 

FIG, 12 is a slide showing a mapping array. 

FIG. 13 is a slide illustrating sinusoidal dewarping. 

FIG. 14 is a slide illustrating correction of nonlinear scan speed. 

FIG. 15 is a slide illustrating the result of a correction. 

FIG. 16 is sUde illustrating divergence correction. 

FIG. 17 is a slide illustrating refraction at interfaces. 

FIG. 18 is a combination of all techniques showing flie images which can be achieved. 
FIG. 19 is a slide illustratmg inserted zoom. 

FIG. 26 is a slide illustrating temporal average and speckle reduction. 

F5Gr-21-4ausfeat^s-aflrever4ay^ehmquer 

FIG. 22 is a flow chart illustrating the overlay technique. 

FIG. 23 illustrates OCT images of the temporal anterior chamber angle of a human eye, 
imaged in vivo at 8 Q>s in different stages of dewarping. 
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FIG. 23a (i) is an image of an Intralipid© drop on a coverslip. Notice the bending of the 
flat slip surface and the bump well below the drop. 

FIG. 23a (ii) illustrates a corrected image with a flat surface and no bump, FIG. 23b is a 
raw image and FIG. 23c illustrates divergence coirection of a handheld scanner and FIG. 24d 
illusti-ates the image of 24b corrected for refraction at the air-cornea boundary (ncomea = 1-38) and 
at the endothelium-acqueous boundary (nacqueous - 1-33). 

FIG. 24 is a sequence of images illustrating the image correction. 
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Detailed Description 
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Real-time imaging and image processing disclosure 



In this sac-i-Torv^ we demonstrate and explain software techniques that are o^se. Jl for 
real-time imaging in OCT (Optical coherence tomography). Special attention is given to 
image processing methods that correct geometric and angular image distortions. 
Examples are geometric distortions due to the applied scanning configuration, nonlinear 
scanner movements and refraction. Polar scanning configurations, useful for endoscopy 
are treated. We are able to practice both . online and ofOLine processing, because we are 
able to save the complete acquired digital data. A special history function greatiy 
simplifies and sometimes enables to save the best images for documentation and later 
processing. 

1 Image transformations, corrections of nonlinearities, 
achieving geometrical and angular correct images 

For display, P' = (x',y') in the acquired data, the raw image, has to be transformed into 
P = (x.y) in the target image. In principle, this can be done in forward (P = f(P')) or 
backward (P' = F(P)) direction. For forward mapping the target position for a given data 
pOTnrisxalcuiated, This±as"iaricey-disativantageHSince^ 

be between target pixels, sophisticated algorithms have to be applied to distribute its 
value onto the neighboring pixels to prevent dark spots and ambiguous assigned pixels, 
which leads to a high computational expense. Backward mapping avoids this 
disadvantage by mapping each target pixel to a location in the acquired image, then using 
simple interpolations to obtain its value. If the backward transformation is fixed, it can be 
implemented with^Iookup table to achieve real-time imaging (Sect 2). In the raw image. 
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x' and y* denote the coordinates across and along A-scans (single depth scans). To obtain 
the brightest possible images with OCT, the field of view with a physical width w and 
depth d is centered a focal length f away from the lens on the optical axis. The size of the 
raw (data) and target (display) image is n' x m' (h x v) and n x m, respectively. The 
backward transformation, separated in both coordinates, is denoted as 

x'=F,(x,y) 
y = F,(x.y) 

!•! Correct aspect ratio 

Often the width and height of acquired data points, given in physical dimensions they 
represent, are not equal. Direct mapping of this raw data acquired on the screen would 
lead to the wrong aspect ratio in display. The following transformations correct this: 

" (2) 

> _ m 

with the condition that the target image has the correct aspect ratio 

m = n-^. (3) 
w 

1.2 Bidirectional scanning 



Fig, I 
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When an OCT system utilizes double-sided scanning (i.e., A-scan acquisition during both 
directions of the reference arm scan), a transformation is necessary to rectify the 
alternate, reversed A-scans (Fig. 1). Again, a static backward transformation can be 
formulated to transform the acquired image array into the image array to be displayed. 

x/2 if round|^y ^ j is even 
-x/2 else (4) 



'ijidirect.x (^» V) ~ 

Fbcdimct.y(x.y) = y/2 



With bidirectional acquisition (see also Sect L3 and 3.1) and a duty cycle t] lower 
100 %, there Mead' pbcels acquired between the end of the forward and the begin of the 
backward scan, which have to-be omitted in Eq, (4). This correction Avili be simplified if 
the duty cycle is center in the physical forward and backward scan. 
Compare also: "Chapter 6. System Integration and Signal/Image Processing", 
Sect. 6.2.2.1.5 

1.3 Registration 

When implementing double-sided scanning correction, it is necessary that consecutive A- 
scans are registered with respect to one another. Depending upon the OCT system 
configuration, the scan registration may be accomplished by adjusting the delay line 
optics. In a clinical setting, however, where the hardware is closed, it may be more 
desirable to impl«nent a software registration mechanism. This could be accomplished 
by allowing manual adjustment by the operator, or by an automatic registration 
algorithm. A sunple way, without computational e;q)ense, is to change the position of the 
start pixel of the acquired scan on the fiamegrabbea: card within the window of a comlete 
linescan. Because this shifts the position of the forward and the backward scan by 1 pixel, 
the registration can only be done with an accuracy of 2 pixels. Fine adjustments result in 
changes in the backward transformation, which can be precalculated in the mapping array 
(Sect. 2). Algorithms for automatic registration will be discussed in Sect. 3.1), 

1.4 Transpose \* > 

Raw images are normally acquired A-scan per A-scan with a framegrabber, forming line 
by line in the raw image. The columns in the raw image represent different depths. For 
display, subsequent A-scan form column by column on the screen, a transpose operation 
is therefore necessary: 
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'"transpose.x (^^y) ^ .(5) 
f^lranspose.y (^» Y ) ~ ^ 

In a similar way, arbitrary rotations are possible. 
1,5 Polar image transformation 

This technique is used in combination with a rotational scanning probe (e.g. an 
endoscopic probe). With this probe A-scans are taken in a radial fashion, with the probe 
constantly rotating. Therefore x* and y' are rather polar coordinates: 

n (6) 
e(y') = 27t-5^ 

with the radius rp«,be of the probe and the imaging depth d. R and e are dimensioniess. 
They can also be expressed in target coordinates 

n/2 d (7) 



e(x.Y) = arctan|^^j 



By combining Eq. (6) and (7) the backward transformations can be obtained 

d/2 d d 



fra<flal,x(^>y)~ 

F«aiau(x.y) = ^«'-ctan[^j 



(8) 



with Wrad the circumference lengdi of the acquisition. 

Compare also: "Chapter 6. System Intep fration and Signal/Image Processing". 

Sect 6,2.2.1.4 and "Multistep image dewarping of OCT images using Fermat's principle 

and moping arrays". Sect. 1.2, 

1.6 Correction of nonlinearities in scanning mirror movements 

An acquired OCT image will be warped if the spatial distribution of the acquired data 
does not directly correspond to the spatial distribution of scattering profile of the sample. 



)_„030117a4A2J.> 



wo 03/011764 



12 



PCT/US02/24y21 



This occurs in OCT imaging when the image data is not sampled at regular intervals in 
space. For example, if the scanning motion of the OCT probe or delay line is not linear 
with time, and the data is sampled at regular intervals in time, then the image will be 
warped. If the scan nonlinearity is a known function of time, however, the image can be 
'de-warped' by an appropriate spatial transformation. This is the case, for example, for the 
sinusoidal motion of a resonant scanning device. In this case, the coordinate 
corresponding to the resonant scanner can be transformed by a sinusoidal function with a 
period corresponding to the period of the scan in image space. Alternatively, if an 
accurate reference signal is available, a corresponding sampling trigger signal could be 
generated to sample nonlinearly in time such that the image is sampled linearly in space. 
This latter technique is common in Fourier transform spectrometers, and has previously 
been applied in high-accuracy interferogram acquisition in OCT. 

In this disclosure we propose a method to dewarp the image accurately after acquisition 
using the position sensor information recorded before OCT imaging. 
This correction of image distortions due to scanner nonlinearities works in both, cross- 
sectional imaging (like in standard OCT) and in en face unaging (like in OCM (optical 
coherence microscopy) or scanning confocal microscopy). 

1.6.1 Fast axis scan 

The fast axis scan is the most likely to show nonlinearities in the scan. When moving or 
turning a mirror at a fast pace effects of momentum and inertia prevent the control 
electronics of the scanner to regulate the scanner position exactly to the sawtooth or 
triangular driving waveform used as command voltage for linear scans. But normally 
scatmers (galvanometers) provide a position sensor output with can be sampled into the 
fiamegrabber input (either instead of the OCT signal or into a different input of the 
fiamegrabber). Gain and offset of the fiamegiabber have to be adjusted to have flie sensor 
input to almost fill but not overfill the fiamegrabber input voltage range. Than for each 
pixel position in the A-scan the corresponding physical position in the sample, given by 
the sensor output, can by recorded. Assuming that the nonlinearities are repetitive 
"Between subsequent scans, averagingTheTpositTorTv^ 

scans reduces the noise by a factor of the square root of the number of scans. The 
standard deviation gives the confidence assigned to this position detection. Let's assume 
an A-scan consists out of m' pixels at the time positions y'i, having position values p'i . 
The p' j are centered around zero, with an amplitude of dbAp. For the backward 
transformation Fscnsor,y(x,y) the algorithm looks for position y'i,min with the minimum 
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p'i.min and the position y'i,max with the maximum p'i,max. With this data set the reverse 
transformation can be calculated from the reverse interpolation (e.g. linear or spUne) 



with p'i limited to (p'umin-p'i^ax) and y'i limited to (y'i,m!n.-y'i,max). 5ce-. . 

B = Interpolation (AA,BB,A) uses the input vectors AA and BB to estimate the value B 

for the position A. 

If the fast axis scanner is used for carrier frequency generation, then strong nonlinearities, 
that can be corrected position wise lead to strong changes in the carrier firequency. 
Therefore either the first bandpass limiting the detectable signal bandwidth has to be 
opened up to pass those wider range of signal frequencies (at the expense of SNR), 
tracking filters or tracking local oscillators have to be used to adapt the current bandpass 
center frequency to the center carrier frequency or a phase modulator for constant center 
frequency have to employed. 

L6.LI Special case: resonant mirror 

For fastest scanning only resonant mirrors can be employed, with scan frequencies of 
several kHz, necessary for real-time OCT imaging. When using a hi^ duty cycle of 
50 % or higher, the scan becomes more and more nonlinear. Because of the operation in 
the resonance mode, the scan follows very exactly a sinusoidal pattern. With fliis 
knowledge the backtransformation can be determined without measuring the position 
waveform. 

The raw image is captured with n' pixels per A-scan and m' A-scans. In principle there 
would be n'pp pixel (peak to peak) available in an A-scan, therefore the duty cycle ii is 
defined as 




The forward transformation from raw'(x\ y') into target (x, y) coordinates is defined as 




(9) 




(11) 
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with the full peak to peak amplitude m'pp of the A-scan in the intermediate coordinates. 
Because r| < 100 % only part of the full sinusoidal motion is visible, the intermediate 
image span m-, pixel in depth, mpp can be calculated from 



2 



rr 

sin 



(12) 



Therefore the back-transformation will be as follows 



Fsinus^(x, y) = ^arcsini ^ j 



(13) 

'"sinus.y Y ) ~ ^ *^ 

See also "Multistep image dewaiping of OCT images using Fermat's principle and 
mapping arrays" Sect, 2. 

1 .6.2 Slow axis scan 

Although less obvious, the scanning nonlinearities in the slow axis scan can be corrected 
similarily to Sect 1.6.1. 

!•? Correction of geometric errors due to scanning configuration 

Depending of where the scanning pivot is located in respect to the final lens (focusing the 
light into the tissue) different scanning regimes can be defined. This different regimes can 
be distinguished into converging, telecentric and diverging scanning configurations. 
Diverging and converging scanning configuration lead to geometric image distortions, 
because the field of view is not rectangular. For this disclosure we propose a 
transf ormatio n tocorrect this distortion. This math is described in "Correction of image 
distortions in OCT based on Fermat*s principle" 
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2 Mapping array 



Fig. 2 

A mapping array has the same dimensions as an expected output image. This array 
represents the point set of an output image in which each element contains the location of 
an input pixel. With the mapping array, the value set of the output image can be obtained 
by backv^'ard mapping to the input image. In a software implementation of image 
tKmsformations, the required mapping array needs to be created only once and stored in 
memory. This approach minimizes computation time while imaging as compared to 
iterative formula-based calculation of image transformations in real-time. Every static 
image transformation can be formulated with this lookup technique, e.g. the correction of 
the aspect ratio (Sect 1.1), bidirectional scanning (Sect. 1.2), registration (Sect, 1.3), 
transposition and arbitrary rotation (Sect 1,4), polar image transformatipn (Sect 1.5), 
correction of nonlinearities in scan mirror movements (Sect 1.6) and the correction of 
geometric errors due to scanning configurations (Sect 1,7), Because all methods are 
based on backward transformations, they can be cascaded, e-g. combining bidirectional 
scamiing with the correction of nonlinearities of scanner movements and conrecdon of 
geometric errors due to the scanning configuration: 

f\otaUy(^« y) ^ Fgeometfic,x.y('^ Mannef,x.y (FbWiiBCl,x,y V))) (^^^ 

with this being a shorthand for the transformations in x and y. 

The mq^ping array consists of pointers ptr(x,y) that point to the direct memory location 
of-the-raw-dat^irhere-given ^er-the-next neighbGr4nterpolatien^ 

ptr{x. y) = Raw_data_start_address + round(Ftote,^(x, y))+ ^^^^ 
round(Ftoto, y (x. y))* bytes _ per _ line 

When data is to be displayed the gray value g(x,y) at the target pixel (x,y) is given by 

g(x.y) = *ptr(x.y) (16) 



OOCID: <WO_0301 17e4A2_L> 



wo 03/011764 



16 



PCT/US02/24'/21 



where **p' denotes the access to the byte vy^here p is pointing to, 
2-1 Inserted zoom 

The mapping array is also capable of real-time zooming of a given window, with no 
penalty in image transformation time. To achieve that, a certain portion (e.g. the upper 
right quadrant,"zoom target") is reserved for the online zoom. Second, is small rectangle 
("zoom source") is defined somewhere else in the image, smaller than the zoom target. 
Now the pointer defined in the zoom target are replaced by pointers that point into the 
raw data for pixels in the zoom source. Since the zoom target in bigger than the zoom 
source, the raw data is sampled finer than in the rest of the image (Fig. 3). Since the zbon 
data is resampled, not just blown up, details that previously were being hidden, become 
visible. This is especially true, when the source image is bigger than the target image, or 
due to strong nonlinearities in the transformation many source pixel are hidden. An 
example for that is the polar transformation, where close to the probe many raw data 
points are not shown (Fig, 3), 



Fig. 3 

3j2-— Interpolation 

When doing the backward transformation, the pixel location for the data requested in the 
raw image will be not an integer position. The fastest method to do the interpolation is 
the "next neighbor" interpolation, which just rounds the source pixel location to the next 
integer position. This method is lossy and generates a blocky target image appearance, 
especially for higher zooms, A better method is using the bilinear interpolation which is 



300CID: <WO__03011764A2_I_> 



wo 03/011764 



17 



PCT/US02/24721 



more computalional expensive. We propose a method to optimize this method to still 
allow real-time imaging. Each entry ptr(x,y) in the mapping array is therefore extended 
by 4 bytes wi to W4 with the relative weight of the neighboring four pixel, coded in the 
range of 0.,255. The target of the ptr(x,y) is always rounded down with the floor() 
operation. 

ptr(x, y) = Raw_data_start_address -h floor(Ftota,^(x, y))-f 
f loor(Ftoia,,y (x. y)) * bytes _ per _ line 



(17) 



When data is to be displayed the gray value g(x,y) at the target pixel (x,y) is given by 

'^*ptr(x,y)*w^+ ^ 
*(ptr(x,y) + 1)*W2 + 

* (ptr(x, y ) + bytes _ per _ line) * W3 + 

* (ptr(x. y) + bytes „ per _ line -h 1) * y 



g(x.y) = 



shr2 (18) 



shr 2 devides by 4 (but is fester) and nonnalizes the range of g back to the original range 
ofthe raw data. 



For real-time color Doppler display we combine our standard mapping array with all its 
advantages with a second lookup table. The data sampled for the structural image is 
bytewise interleaved with the flow data, forming I6bit values of data containing all 
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information for a data point. To dispLay a pixel as a flow value, two conditions have to be 
met: (1) the flow must exceed a tlireshold and (2) the value of the structural image at this 
point has to be high enough to validate the presents of blood, which is highly scattering at 
1 300 mm. Fig. 4 shows the flow chart for this decision process. Checking these two 
conditions while displaying would be too computational expensive for real-time display. 
But by precalculating a lookup table for the resulting color for all possible combinations 
of structural and flow value, this decision process is reduced to a simple memory access 
by using the 16 bit value of the combined structural and flow value as a pointer into this 
lookup table. 

3 Online image processing 

We iinplemented a couple of other important online processing methods to improve 
image quality. 

3.1 Auto registration 

When using bidirectional scanning, the forward and backward scan have to be carefully 
registered to achieve a sharp image. The easiest way to implement the registration is by 
user input. Unfortunately, the registration has to be readjusted about every 5- 10' minutes • 
of imaging, due to small instabilities in the linesync generation. Manual registration is 
been complicated by the fact that the human observer only can detect that the image gets 
*fuzzy% but not clear in which direction to change the registration without trying. 
Secondly, in a moving sample it is difficult to correctly see any small misregistration. An 
example is endoscopic imaging: In air, the outer sheath can act as a reference line for 
registration, but inside the patient this reflection vanished due to tissue contact and 
variable amounts of fluid on the probe. 

To address this problem, we propose utilizing the reflection at an unaffected surface like 
the inner sheath surface of the endoscopic probe or if a glass plate is present in the 
sample arm, using the back surface of this glass plate. Since this reference surface is 
always ihe-sam^distancc^vvay-^-w^-can^^ 

backward scan, limited to the depth range, where the reference line is expected. The 
computer than calculates the cross-correlation function for a full range of different lag 
and changes the registration to achieve maximum cross-correlation at zero lag. To 
prevent oscillations, this is done only every few images with following algorithm: If the 
calculated maximum lag is big, the next registration step will reduce to the half. If the 
cross-correlation has its maximum close to zero lag, the registration will be only change 
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by the smallest possible step. Experiments showed that within seconds the registration 
reached its best point and stayed there. 

3.2 Running average with lateral regisf ration 

To improve the SNR and to reveal details deeper into the tissue we applied an averaging 
teclmique to subsequent images. This was implemented as a running average to sustain 
real-time imaging 

Sn = ^n-i 'n ~ 'i-a /-jqn 

D„=S„/a 

with the n*^ sum Sn, the n*** image acquired I„, the n^ displayed image D„, and the running 
average length a. Lateral motion artifacts were strongly reduced by lateral registration 
through cross-correlation of lateral stripes of 20 pixel in depth in a imaging where sample 
signal is to be expected. 

3.3 Online detection of front surface 

For depth gain compensation and correction of the index of refiraction, as well as 
correction of refiraction, it is useful to online detect the topmost surface of the sample, 
because there the absorption and index of refraction changes dramatically at this 
boundary. We accomplished this detection by using following steps 

• Gaussian filter of with a larger axial size thlan lateral size. 

• Edge detection with a axial Sobel jSlter with the length of the axial Gaussian jSlter, 

• look for liie maximum in the edge image axial every 30 A-scans. 

• Use this points as starting points for active contours. 

3.3.1 Correction of index of refraction below surface with dynamic 
mapping array 

Since OCT measures the optical pathlength rather than the physical distance, OCT 
images taken in air are distorted in a way that the tissue layers appear ttiicker than they 
are. This error is most obvious with samples that have wavy surfaces, like skin with 
demial ridges. These waves will have artifacts in deeper layer structures, disturbing the 
diagnosis (Refraction is another issue that will be discussed in Sect 5.1). We purposes 
here the application of a mapping array with a dynamic added ID overlay. This overlay is 
chosen depending on the axial position in which the surface v/as detected in this A-scan 
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and the index of refraction and allows real-time display with correction of index of 
refraction. 

3-3.2 Depth gain compensation 

The strength of the OCT signal degrades exponentially with depth. Since usually the 
OCT signal is displayed on an exponential scale on the screen, this means, the intensity in 
gray values drops linearly from the depth of the scan where it hits the surface. It the top 
surface is known, a linearly with depth growing offset can be added to the OCT signal to 
compensate for the loss. This is limited by the amplification of noise outside the tissue 

4 Streaming, online documentation 

The standard procedure before streaming was to have a live display, and when an image 
with relevant content appears, to jfreeze the display and save this image. Alternative the 
screen output is recorded to S-VHS. Both approaches have strong disadvantages. When 
the feature of interest was only transient on the screen, freezing was maybe to slow to 
keep the image, it is then lost. Recording on VCR reduces the image quality dramatically. 
We propose to record the complete raw digital data to harddrive while imaging (as a 
stream of data). This approach has several advantages 

• retaining the fiill, complete raw data allow for offline processing the data which 
algorithms not available in real-time. 

• Zooming in post-processing shows data not visible during live imaging (c,f. Sect 2. 1) 

• Tlie frames of a stream are appended to an expanding multiple TIFF-file. TIFF is a 
very flexible structure, which also allows for storage of additional information with the 
images, like patient and study information, acquisition parameters, labels and tissue 
identifications and classifications 

When saving as a multiple TIFF, streaming of all acquired images into a circular buffer 
(called the history buffer) for the images worth 10-20 s. We propose this circular buffer 
as a very useful tool to retrieve good images and low load on the hard drive capacity. 
After freezing the acquisition the user has access to all images of tiie last 1 0 to 20 sec 
with hotkeys or by mouse selection. Functions available are going framewise forward or 
backward in the history, cyclic playing of the history images. There is a function to save 
the currently displayed frame or all frame from this history buffer. Hotkeys can be 
associated with VCR like keys as easy mnemonics. Before saving images can be 
classified, the visible organ with shortcut buttons or typing specified, features visible can 
be labeled onscreen with an overlaying label (nondestructively). All this extra 
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information will be saved within single TIFF's, Altemativeiy, all single images, save 
history buffer, and direct streaming will be saved in one file. The idea is to have one file 
per procedure or patient, for easy documentation. All images have a timestamp with a 
resolution of ms saved with it for easy and unique identification. 

5 Offline image processing 

Since we save the complete raw digital data acquired, we can do offline processing 
without any losses, with the advantage of being able to do computational intensive 
calculations to improve image qusdity and feature visibility 

5.1 Correction of index of refraction 

cf. Correction of image distortions in OCT based on Fermat's principle 

5.2 Offline registration 

Even though automatic registration of forward and backward scans during processing 
works reasonable well, best results are obtainable by human intervention in post 
processing, aligning it to a half a pixel precision 

53 Stream view 

To handle the huge amount of images 
1 line/frame 

Info 1 and 2 on both sides and short bars in the middle for fast orientation 
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System Integration and Signai/lmage Processing 



This s«^ianr addresses the integration of the component technologies of OCT, as described 
in the preceeding 5e£:^'»<^*ts., into a complete imaging system (see figure 6,1). This includes both 
hardware considerations including optimal interferometer topologies and scan synchronization 
dynamics, as well as software issues including image acquisition, transformation, display, and 
enhancement. A limited discussion of first-generation slow-scan systems (<1 image/sec) is 
included where it is illustrative; however, most of the discussion centers upon state-of-the-art 
OCT systems acquiring images in near real time. 
6.1 Hardware Implementation 

6.1.1 Interferometer topologies for optimal SNR 
The original and most common interferometer topology tised in OCT systems is a simple 
Michelson interferometer, as depicted in earlier chapters. In this design, low-coherence source 
light is split by a 50/50 beamsplitter into sample and reference paths. A retroreflecting variable 
optical delay line (ODL) comprises flie reference arm, while the sample specimen togetibier with 
coupling and/or steering optics comprise the sample arm. Light retroreflected by the reference 
ODL and by the sample is recombined at the beamsplitter and half is collected by a photodetector 
in the detection arm of the interferometer. Heterodyne detection of the interfering light is 
achieved by Doppler shifting the reference light with a constant-velocity scanning ODL, or by 
phase modulating either the sample or reference arm. Single-mode fiber implementation of the 
interferometer has the advantages of simplicity and automatic assurance of the mutual spatial 
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coherence of the sample and reference light incident on the detector. Although this design is 
intuitive and simple to implement, it is apparent that due to the reciprocal nature of the 
beamsplitter, half of the light backscattered Ifrom the sample and from the reference ODL is 
returned towards the source. Light returned to the source is both lost to detection and may also 
compromise the mode stability of the source. Further, a detailed analysis of the signal-to-noise 
ratio in low-coherence interferometry [1, 2] mandates that in order to optimally approach the 
shot-noise detection limit for the Michelson topology, the reference arm light must be attenuated 
by several orders of magnitude (depending ui>on the source power level and the detector noise 
figure). Thus, in typical OCT implementations, up to 75% (or 6 dB) of the optical power 
suppUed by the source does not contribute to image formation. Since light sources with high- 
spatial coherence and low temporal coherence suitable for high-speed OCT imaging in very low 
backscattering samples such as biological tissue are very expensive, optical power is a 
commodity well worfh conserving in OCT systems. 

Using optical circulators, unbalanced fiber couplers, and balanced heterodyne detection, 
two authors have recently demonstrated a new family of power-eflBcient fiber-optic 
interferometer topologies which recover most or all of the light lost in the Michelson 
configuration [2, 3]. The three-port optical circulator is a non-reciprocal device which couples 
light incident on port I to port 11, and light incident on port 11 to port in. Current commercial 
fiber-coupled devices specify insertion losses less than 0.7 dB (I to II, II to III) and isolation (HI 
to II, II to I) and directivity (I to III) greater than 50 dB. Single mode wideband fiberoptic 
couplers are commercially available with arbitrary (unbalanced) splitting ratios. Balanced 
heterodyne reception is an established technology in coherent optical communications [4, 5], and 
has previously been used in OCDR [6] and in OCT [3, 7-9]. 
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Tliree types of new interferometer designs described by Rollins et al [2] utilizing these 
enabling technologies are illustrated in figure 6.2. The first design (figure 2A) uses a Mach- 
Zelinder interferometer with the sample located in one arm and the reference ODL in the other 
arm. The first coupler is unbalanced with a splitting ratio chosen to optimize SNR by directing 
most of the source light to the sample. Light is coupled to the sample through an optical 
ckculator so that the backscattered signal is collected by the delivery fiber and redirected to the 
second coupler. The reference ODL may be transmissive, or alternatively, a retroreflecting ODL 
may be used with a second circulator. Design Aii of figure 6.2 is similar to Ai except that instead 
of balanced heterodyne detection, a second unbalanced splitter is used to direct most of the 
sample light to a single detector. Although less expensive to implement, the performance of the 
single detector version is significantly worse than the balanced detector version since a single 
detector does not suppress excess photon noisc; 

Interferometer design B is similar to design A, as shown in the schematics labeled Bi and 
Bii in figure 6.2. In this case, a retroreflecting ODL is used without the need for a second optical 
circulator. Configuration Bii has recently been demonstrated for endoscopic OCT [3]. 

Design C uses a Michelson interferometer e£Bciently by introducing an optical circulator 
into the source arm instead of the sample arm, as in designs A and B. Configuration Ci utilizes a 
balanced receiver. Configuration Ci has the significant advantage that an existing fiber-optic 
Michelson interferometer OCT system can be easily retrofitted with a circulator in the source 
arm and a balanced receiver with no need to disturb the rest of the system. One drawback of 
configuration Ci is that more light is incident on the detectors than in the other configurations. 
Although the balanced receiver is effective in suppressing excess photon noise, a lower gain 
receiver is necessary to avoid saturation of the detectors. In a high speed OCT system, however. 
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this is not an issue because a lower gain receiver is necessary to accommodate the broad signal 
bandwidth. Design Ci has also recently been demonstrated for use in endoscopic OCT [9], 
Design Cii uses an unbalanced splitter and a single detector. 

A published theoretical analysis of the signal-to-noise ratio for all of the designs illustrated 
in figure 6,2 [2] illustrates that with the proper choice of coupler splitting ratios and balanced 
receiver configurations, the three dual-balanced detection configurations may provide an increase 
in SNR of between 6 and 8 dB as compared to the standard Michelson topology. The expected 
increase in SNR beyond the 6 dB power loss of the Michelson is due to ttie additional capability 
of balanced reception to reduce excess photon noise [6]. The single-detector versions provide a 
more modest gain in SNR. 

6.1.2 Analog Signal Processing 

The critical function of the analog signal processing electronics in an OCT system is to 
extract the intcrfeiometric component of the voltage or current signal provided by the detection 
electronics with higji dynamic range and to prepare it for analog-to-digital conversion. Other 
functions which could potentially be performed at this stage include signal processing operations 
for image enhancement, such as deconvolution, phase contrast, polarization-sensitive imaging, or 
Doppler OCT, Although to date almost all such image enhancement processing has been 
perforaied in software, as high-speed imaging systems become more prevalent, pre-digitization 
time-domain processing will inevitably become more sophisticated. For example, a recent 
demonstration of real-time signal processing for Doppler OCT have utilizjed an analog approach 
[10]. 

Demodulation of tiie interferometric component of the detected signal may be perforaied 
using either an incoherent (i.e. peak detector) or coherent (quadrature demodulation) approach. 
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Many early low-speed OCT systems for which the Doppler shift frequency did not exceed 
several hundred kHz utilized a simple custom circuit employing a commercially available 
integrated-circuit RMS detector in conjunction witii a one- or two-pole bandpass filter for 
incoherent demodulation [11, 12]. Even more simply, satisfactory coherent demodulation can be 
accomplished using a commercial dual-phase lock-in amplifier without any other external 
components [13, 14]. By supplying the lock-in amplifier with a reference sine wave at the 
calculated reference arm Doppler shift frequency and using the lock-in time constant as a good- 
quality high-order low-pass filter, the amplitude of the sum of the squares of the in-phase and 
quadrature outputs provides a high-dynamic range monitor of the interferometric signal power. 
In more recent high-speed OCT system implementations employing modulation fi:equcncies in 
the MHz range, more sophisticated electronics based on components designed for the ultrasound 
and cellular radio communications markets have been employed [3, 9, 15-18]. 
6.1.2.1 Dynamic Range Compression 
As discussed in previous chapters, typical OCT systems routinely achieve sensitivity 
(defined as the minimum detectable optical power reflectivity of the sample) well in excess of 
100 dB. Since the electronic interferometric signal amplitude (current or voltage) is proportional 
to the product of the reference and sample arm electric fields and thus to the square root of the 
sample power reflectivity, the signal-to-noise ratio of the electronic signal amplitude usually 
exceeds 50 dB. This value both exceeds the dynamic range of the human visual system (which 
can sense brightness variations of only about 3 decades in a particular scene) and also approaches 
the dynamic range limit of many of the hardware components comprising the signal 
detection/processing/digitizatipn chain. For example, the dynamic range of an A/D converter is 
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given by 2^^ (= 6N dB), where N is the number of bits of conversion; thus an 8-bit converter has 
a dynamic range of only 48 dB. For early OCT systems employing slow pixel digitization rates 
of only a few hundred kilohertz, this latter issue was not a limiting factor since high dynamic 
range A/D converters (i.e. up to 16 bit) are common and inexpensive at these data rates. For 
high-speed OCT imaging, however, MHz digitization rates are required and the dynamic range of 
the digitization step becomes an increasingly important factor both in terms of digitizer cost as 
well as in downstream computation speed. 

In order for an OCT image to be rapidly digitized and meaningfully observed, a means of 
hardware or software dynamic range compression is often employed. This is accomplished by 
transforming the detected sample reflectivity values with a nonlinear operation that has a^ 
maximum slope for low reflectivity values and a decreasing slope for increasing reflectivity - 
values. The obvious and convenieiit method is to display the logarithm of reflectivity in units of 
decibels. The logarithm operation demonstrates the desired transform characteristic, and decibels 
are a meaningful, recognizable unit for reflectivity. The logarithm is not the only possible 
dynamic range compression transform. For example, the ji-law transform of communications 
systems, or a sinusoidal transform could be used, but up to the present time, logarithmic 
compression is universal in display of OCT images. 

Since the first reports of OCT, images have conventionally been displayed in logarithmic 
format, usually having been transformed in software [19]. Although plots of A-scan data can be 
well visualized on a linear scale, in two-dimensional images displayed with a linear intensity 
scale only the brightest features are perceived (see, for example, [20]). It is notable that the 
common definition of axial resolution in OCT as being given by half the coherence length of the 
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light source is valid only when the data is presented on a linear scale; the 3-dB point of a sharp 
reflection in a logarithmically transformed image depends upon the dynamic range of the image 
data and is thus not well defined, but is clearly wider than the half-coherence length. 

Several options for hardware analog dynamic range compression are available. An 
approach which has been utilized in several OCT implementations compresses the dynamic 
range of the signal before sampling by using an amplifier with a nonlinear gain characteristic [9, 
11]. In this way, commonly available data acquisition boards and firame grabbers with linear 
quantum levels can still be used for digitization. Demodulating logarithmic amplifiers are 
commercially available which compress the dynamic range of the OCT signal and also perform 
envelope detection. Alternatively, A/D converters are available with non-linearly spaced 
quantum levels, for sample foUowiQg tiie n-law transform. TTiis would allow the low- 
reflectivity range to be sampled with high A/t) dynamic range, while- sacrificing A/D dynamic 
range vsrhen sampling the high-reflectivity range, where it is not needed. Devices are also 
available which perform sinusoidal transformations. 
6.1.3 Synchronization electronics 
As illustrated in figure 6.1, every OCT system implementation includes at least some fonn 
of optical delay line, sample scanning optics, and digitization/display electronics whose dynamic 
fimctions must be coordinated to acquire meaningfiil image data. In addition to the 
synchronization of these elements which is common to all OCT systems, specially designed 
systems may also require coordination of dynamic properties of flie optical source (e.g., 
frequency-tunable source implementations [21]), detection electronics, or analog signal 
processing electronics (e.g., frequency-tracking demodulators [22]). 
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A diagram illustrating the timing relationships between the elements in a standard minimal 
system is illustrated in figure 6.3, For most common systems which employ depth-priority 
scanning (i.e., in which depth A-scans are acquired more rapidly than lateral scans), individual 
image pixels are acquired in the IkHz — lOMHz range, the reference optical delay has a 
repetition rate in the 10 Hz -lOkHz range, and the lateral sample scanning optics repeat at 0.1-10 
Hz frequencies. The optical delay line is driven by a waveform which is optimally a triangle or 
sawtooth to maximize the duration of the linear portion of the scan and thus the usable scan duty 
cycle, although harmonic and other nonlinear delay waveforms have been used for the fastest 
delay lines yet reported [18, 23]. In either case, the synchronization electronics provide a frame 
sync signal synchronized to the sample lateral scan to signal to the image acquisition electronics 
to start image frame acquisition. Similarly, the synchronization electronics provide a line sync 
signal synchronized to the depth scan to signal the image acquisition electronics to start A-scan 
digitization. In most OCT system implementations described to date, a pixel clock is generated 
by a synthesized source (i.e. by a function generator or on-board A/D conversion timer) at a 
digitization rate given by the line scan rate multiplied by the number of desired pixels per line. 
However, an OCT system specially designed for coherent signal processing (utilized both for 
OCT image deconvolution [24] and for spectroscopic OCT [25]) has been reported which 
utilized a helium-neon based reference arm calibration interferometer to provide a pixel clock 
sync signal coordinated to the actual reference optical delay with nanometer accuracy. Lateral- 
priority scarming OCT systems have also been reported; in this case the timing of the depth and 
lateral scans are reversed [13, 26-28]. 

The hardware comprising the synchronization and image acquisition electronics may be as 
simple as a multifunction data acquisition board (analog-to-digital, digital-to-analog, plus timer) 
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residing in a personal computer. Alternatively, as described in the next section, a standard video 
frame grabber board may be programmed to perform the same functions at much higher frame 
rates. 

6.1.4 Example of an Integrated OCT System 
As an example of an integrated OCT system, a block diagram of a rapid-scan system 
designed for endoscopic evaluation of early cancer is provided in Figure 6.4 [9]. The high speed 
OCT interferometer is based on a published design [18]. It includes a high-power (22 mW), L3 
pm center wavelength, broadband (67 nm FWHM) semiconductor amplifier-based light source, 
and a Fourier-domain rapid-scan optical delay line based on a resonant optical scanner operating 
at 2 kHz. Both forward and reverse scans of the optical delay line are used, resulting in an A- 
scan acquisition rate of 4 kHz. Image data is digitized during the crater two-thirds of the 
forward and reverse scans, for an overall scanning duty cycle of 67%. 

In this system, OCT probe light is delivered to the r^ion of interest in the lumen of the GI 
tract via catheter probes which are passed through the accessory chaimel of a standard GI 
endoscope. A specialized shafiE, which is axially flexible and torsionally rigid, mechanically 
supports the optical elements of the probe. The probe beam is scanned in a radial direction nearly 
perpendicular to the probe axis at 6.7 revolutions per second (the standard frame rate in 
commercial endoscopic ultrasound systems) or 4 revolutions per second. The converging beam 
exiting the probe is focused to a minimum spot of approximately 25 p.m. 

Optical signals returning from the sample and reference arms of the interferometer are 
delivered via the non-reciprocal interferometer topology (figure 6.2Ci) to a dual-balanced 
InGaAs differential photoreceiver. The photoreceiver voltage is demodulated and dynamic range 
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compressed using a demodulating logarithmic amplifier. The resulting signal is digitized using a 
conventional variable scan frame grabber residing in a Pentium II PC. The line sync signal for 
the frame grabber is provided by the resonant scanner controller, the frame sync signal is derived 
from the catheter probe rotary drive controller (1 sync signal per rotation), and the pixel clock is 
generated internally in the frame grabber. 

The PC-based EOCT imaging system is wholly contained in a single, mobile rack 
appropriate for use in the endoscopic procedure suite. The system is electrically isolated and the 
optical source is under interlock control of the probe control unit The system meets institutional 
and federal electrical safety and laser safety regulations. The data capture and display subsystem 
acquires image data at a rate of 4000 lines per second using the variable scan franie grabber. 
Alternate scan reversal is performed in software in order to utilize both forward and reverse scans 
of the optical delay line, followed by rectangular-to-polar scan conversion using nearest-neighbor 
interpolation (see below). Six hundred (or 1000) A-scans are used to form each image. A 
software algorithm performs these spatial transformations in real time to create a full-screen 
(600x600 pixels) radial OCT image updated at 6.7 (or 4) frames per second. Endoscopic OCT 
images are displayed on the computer monitor as well as archived to S-VHS video tape. Foot 
pedals controlling freeze-fi:ame and frame capture commands are provided, allowing the 
endoscopist to quickly and effectively acquire data using the system. 
6.2 Image Acquisition and Display 

6.2.1 High speed data capture and display 
6.2.1.1 Franie grabber tectmology 

Frame grabbers are designed to digitize video signals, such as from CCD cameras, CID 
cameras, and vidicon cameras. If each frame of video signals is 640x480 pixels, the amount of 
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memory needed to store it is about one quarter of a megabyte for a monochrome image having 8 
bits/pixeL Color images requires even more memory, approximately three times this amount. 
Without a frame grabber, most inexpensive general-purpose personal computers cannot handle 
the bandwidth necessary to transfer, process, and display this much information, especially at the 
video rate of 30 frames per second. As a result, a frame grabber is always needed in an imaging 
system when displaying images at or approaching video rate. 

A block diagram of a simple frame grabber is shown in Figure 6.5. Typical frame grabbers 
are functionally comprised of four sections: an A/D converts, programmable pixel clock, 
acquisition/window control vinit, and frame buffer. Video input is digitized by the A/D converter 
with characteristics, such as filtering, reference and offset voltages, gain, sampling rate and the 
source of sync signals controlled programmatically by the programmable pixel clock and 
acquisition/window control circuits. The frequency of the programmable pixel clock determines 
the video input signal digitization rate or sampling rate. In addition, the acquisition/window 
control circuitry also controls the region of interest (ROI) whose values are determined by the 
user. Image data outside of the ROI is not transfoxed to the frame buffer, and not displayed on 
the screen. 

In order to utilize a frame grabber for general-purpose high-speed data acquisition, it is 
important to understand the nature of video signals and how their acquisition is controlled by 
frame grabber driver software. A video signal comprises a sequence of different images, each of 
which is refiared to as a frame. Each frame can be constructed from either one (non-interlaced) or 
two (interlaced) fields, depending upon the source of the signal. Most CCD cameras generate 
interlaced frames* The even field in an interlaced frame would contain lines 0, 2, 4, the odd 
field would contain lines 1, 3, S, and so on. 
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Figure 6.6 illustrates the components of a single horizontal line of non-interlaced video as 
well as the visual relationship between the signal components and the setting of the 
corresponding input controls. Upon the same basis. Figure 6.7 shows the components of a single 
vertical field of video as well as the relationship between the signal and the setting of the 
corresponding input controls. 

6.2.1.2 False color/gray scale mapping 

Once digitiz:ed by a conventional A/D converter or frame grabber, two-dimensional OCT 
image data representing cross-sectional or en face sample sections is typically represented as an 
intensity plot using gray-scale or false-color mapping. The intensity plot typically encodes the 
logarithm of the detected signal amplitude as a gray scale value or color which is plotted as a. 
function of the two spatial dimensions. The choice of the color mapping used to represent OCT 
images has a significant effect on the perceived impact of the images and on the ease (and 
expense) with which images can be reproduced and displayed. 

Many authors [13, 14, 18, 21, .29-31] have used the standard linear gray scale for OCT 
image representation, with low signal amplitudes mapping to black and strong reflections 
appearing as white. Some groups [3, 32, 33] have opted for a reverse gray scale with strong 
reflections appearing as black on a white background, primarily motivated by the relative ease of 
printing and reproducing such images as compared to the standard grsiy scale. A large variety of 
false color maps have been applied to OCT images; the most widely used being the original blue- 
green-yellow-red-white "retina" color scale designed by David Huang while a graduate student at 
MIT for inclusion in the original paper describing OCT in 1991 [19]. This color was 
adopted by Humphrey Systems, Inc. as the exclusive map for their OCT retinal scanner, and has 
also been used in some presentations of non-ophthaimic data [19]. A plot of the retinal felse- 
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color scale in RGB color space is reproduced in figure 6.8, and a comparison of an in vivo 
liuman retinal OCT image in each of the three color scales is provided in figure 6.9. 

Recent innovations in OCT technology which enhance structural or morphological imaging 
by the addition of functional information (e.g. Doppler flow imaging [34, 35], polarization- 
scnsiiive OCT [36], spectroscopic OCT [25, 37J) face a difficult problem of representing all of 
the resulting information in a single image. To date, two different pseudo-color image coding 
schemes have been employed to combine- depth-resolved blood flow imaging or spectroscopy 
with conventional OCT reflectivity imaging. These are based on two standard color models 
described in image processing texts [38]. Briefly, these are the RGB (red, green, blue) and HSL 
(hue, saturation, luminance) models. Hue is associated with the perceived dominant wavelength 
of the color, saturation is its spectral purity, or the extent to which the color deviates from white, 
and luminance is the intensity of color. In the RGB model, the relative contributions from red, 
green, and blue are used to describe these properties for an arbitrary color. In the HSL model, 
color intensity is controlled independently from the hue and saturation of the color. 

Doppler OCT imaging [34, 35] has adapted an RGB color map to simultaneously indicate 
reflectivity and Doppler shifts. Blood flow data are thresholded to remove noise and 
superimposed on the reflectivity image. The standard linear gray scale is used to represent 
amplitude backscatter, v^^ereas blood flow direction is indicated with hue (red or blue for 
positive or negative Doppler shifts, respectively). Higher luminance indicates increased flow 
magnitude. 

Recently, a variation of the HSL colormap was introduced [37] for combining backscatter 
intensity and depth-resolved spectral shifts in tissue. As with the previous scheme, hue denotes a 
shift in the backscatter spectrum, where red, green, and yellow designate positive, negative, and 
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no spectral shift, respectively. Saturation of each hue indicates tissue reflectivity, and the image 
contains constant luminance. 

6.2.2 Image transformations 

In general, images may be defined in terms of two elementary sets: a value set and a point 
set [39]. The value set is the set of values which the iniage data can assume. It can be a set of 
integers, real, or complex numbers. The px^int set is a topological space, a sub-set of n- 
dimensional Euclidean space which describes the spatieJ location to which each of the values in 
the point . set are assigned. 

Given a point set and a value set F, an image / can be represented in the form 
/ = {(Jc,a(x)) : X e JSr,a(jc) e F} . (6,1) 
An elem^t of the image, (x, a(x)), is called a pfacel, x is called the pixel location, and a(pc) is flie;; 
pixel value at the location X. 

Two types of image transformations are of interest in processing of OCT images. Spatial : 
transformations operate on the image point set and can accomplish such operations as zooming,: 
de-warping, and rectangular-to-polar conversion of images. Value transformations operate on the 
value set and thus modify pixel values rather than pixel locations. Examples of useful value 
transformations include modifying image brightness or contrast, exponential attenuation 
correction, or image de-speckling. 

6.2.2.1 Spatial Transformations 

A spatial transformation defines a geometric relationship between each point in an isiput 
point set before transformation and the corresponding point in the output point set A forvyard 
mapping function is used to map the input onto the output. A backward m£43ping function is used 
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to map tlie output back onto the input (see figure 6.10). Assuming that [u, v] and [x, y] refer to 
the coordinates of the input and output pixels, the relationship between them may be represented 
as 

v]-[A>,vXr(z/,v)] or [u.v] = [Ulix,ylV(x,y)l (6.2) 
where X and Y are the forward mapping functions and U and V are the backward mapping 
functions. 

6.2.2.1.1 Spatial Transformation Matrices 

In a rectangiilar coordinate system, linear spatial transformations (e.g. translation, rotation, 
scaling, and shearing etc.) can be expressed in matrix form using a transformation matrix T, [39] 
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The 3x3 transformation matrix Tj can be best understood by partitioning it into 4 separate 
submatrices. The 2x2 submatrix ^" specifies a linear transformation for scaling, 

shearing, and rotation. The 2x1 submatrix [a,3 ^23]^ produces translation. The 1x2 submatrix 
[flj, J produces perspective transformation. The final 1x1 submatrix [^33] is responsible for 
overall scaling and usually takes a imity value. Note that the superscript T denotes matrix 
transposition, whereby rows and columns are interchanged. Examples of simple useful 
transformation matrices in rectangular and polar coordinates are provided in Table 6.1. 
6.2.2.1.2 Mapping Arrays 
Non-linear spatial transformations which cannot be performed using transformation 
matrices (e.g. coordinate system conversions) can be performed using a mapping array. 
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Specification of a mapping array is also a useful and necessary step in the computer 
implementation of linear spatial image transforms. A mapping array has the same dimensions as 
an expected output image. This array represents the point set of an output image in which each 
element contains the location of an input pixel. .With the mapping array, the value set of the 
output image can be obtained by backward mapping to the input image. In a software 
implementation of image transformations, the required mapping array needs to be created only 
once and stored in memory. This approach minimizes computation time while imaging as 
compared to iterative formula-based calculation of image transformations in real time. 
6.2.2.1.3 Image Rotation 
Image rotation is a commonly used image transformation in high-speed OCT systems in 
which depth-priority OCT images (those acquired using a rapid z-scan and slower lateral scan) 
are captured using a frame grabber (which expects video images having a rapid lateral scan). The 
image rotation transformation is illustrated in figure 6.1 1 , A 90-degree-rotation mapping array is 
created to reconstruct an image of the sample from the frame buffer of the frame grabber. 
In this case, the spatial transformation used for creating a mapping array is 
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where 

/ = (0, 1 j =ra a j), 

m =(0. I n ^(0. 1, ....N),I^N,andJ = N, 

6.2.2.1.4 Rectangular*-to-polar conversion 
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Rectangular to polar conversion is necessary when image data is obtained using a radially 
scanning OCT probe, such as an endoscopic catheter probe [9, 40]. The A-scans will be 
recorded, by a frame grabber for example, sequentially into a rectangular array, but must be 
displayed in a radial format corresponding to the geometry of the scanning probe, as illustrated in 
figure 6.12. The forward mapping operations are: 



where x and y are the rectangular (Cartesian) coordinates and r and G are the polar coordinates. 
The inverse mapping operations are: 



In the calculation of 6, an additional conditional is necessary because of ambiguity between the 
first and titiird, and second and fourth quadrants of the polar coordinate system. 



When an OCT system utilizes double-sided scanning (i.e., A-scan acquisition during bolii 
directions of the reference arm scan), a transformation is necessary to rectify the alternate, 
reversed A*scans (figure 6.13). Again, a mapping array can be constructed to transform the 
acquired image array into the image array to be displayed. When implementing double-sided 
scanning correction, it is necessary that consecutive A~scans are registered with respect to one 
another. Depending upon the OCT system configuration, the scan registration may be 
accomplished by adjusting the delay line optics. In a clinical setting, however, where the 
hardware is closed, it may be more desirable to implement a software registration mechanism. 



x(r,e) = rcos(G) 
y(r,e) = rsin(G)' 



(6.5) 



Kx,y) = Vx^+y^ 




(6.6) 



6.2.2.1.5 Double-sided scan correction 
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This could be accomplished by allowing manual adjustment by the operator, or by an automatic 
registration algorithm. 

6.2.2.1.6 Image de-warping 

An acquired OCT image will be warped if the spatial distribution of the acquired data does 
not directly correspond to the spatial distribution of scattering profile of the sample. This occurs 
in OCT imaging when the image data is not sampled at regular intervals in space. For example, if 
the scarming motion of the OCT probe or delay line is not linear with time, and the data is 
sampled at regular intervals in time, then the image will be warped. If the scan nonlinearity is a 
known function of time, however, the image can be 'de-warped* by an appropriate spatial 
transformation. This is the case, for example, for the sinusoidal motion of a resonant scamiing 
device. In this case, the coordinate corresponding to the resonant scaimer can be transformed by 
a sinusoidal function with a period corresponding to the period of the scan in image space. 
Alternatively, if an accurate ref<^nce signal is available, a correspondiag sampling trigger signal 
could be generated to sample nonlinearly in time such that the image is sampled linearly in 
space. This latter technique is common in Fourier transform spectrometers, and has previously 
been applied in high-accuracy interferogram acqtiisition in OCT [24]. 

6.2.2.1.7 Motion Artifact Reduction by A-Scan Registration 

In OCT imaging of living human or animal subjects, motion arti&ct is a serious concern 
because the inherent spatial resolution of OCT imagmg often exceeds the extent to which the 
subject is capable of remaining motionless imder voluntary control. The most straightforward 
approach to eliminate this problem is to acquire images more rapidly, and the recent advent of 
high-speed scanning systems has in large part alleviated this concern [18, 33]. However, in cases 
where the available signal-to-noise ratio precludes high-speed imaging (as, for example, in 
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retinal imaging where high-power sources cannot be used), image processing techniques can be 
applied. 

The technique of cross-correlation scan registration for retinal image motion artifact 
reduction was developed by Hee and Izatt in 1993 [12, 41]. In this technique, an estimate of the 
patient's axial eye motion during image acquisition was formed from the peak of the cross- 
coirelation Ri(k) of each A-scan in the image with respect to its neighbor: 

Rim^EFiOlFMU-k). (6.7) 

Here, Fi(/) is the longitunidal scan data at the transverse scan index /, where j is the longitunidal 
pixel index. The profile of the motion estimated from the locus of the peaks of Ri(k) was then 
low-pass filtered in order to separate motion artifacts (which were assumed to occur at relatively 
high spatial frequency) from real variations in the patirat's retinal profile (which were assumed 
to occur at relatively low spatial frequency). The smoothed profile was then subtracted from the 
original profile to generate an array of offset values which were applied to correct the positions 
of each A-scan in the image. An illustration of ttds procedure and its results is provided in figure 
6.14. 

It is worth noting that the position of the peak of the cross-correlation fimction in retinal 
images appears to depend heavily upon the position of the retinal pigment epithelium (RPE). In 
images in which the strong RPE reflection is absent in some A-scahs (for example, underneath 
strongly attenuating blood vessels), the cross-correlation technique is subject to registration 
errors. In such cases, a motion profile may alternatively be obtained by thresholding the A-scan 
data to locate the position of a strong reflectivity transition within the tissue structure, such as 
occurs at the inner limiting membrane. Thresholding at this boundary has recently been applied 
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for A-scan registration of Doppler OCT images in the human retina [42]. In this case, the 
velocity data was also corrected by estimating the velocity of the patient motion from the spatial 
derivative of the scan-to-scan motion estimate and from knowledge of the A-scan acquisition 
time. 

6.2.2.1 Value Set Operations 

In contrast to spatial transformations, value set operations modify pixel values rather than 
pixel locations. Spatial filtering using convolution kernels are fundamental to image processing 
and can of course be ^plied to OCT images. Examples of useful convolution kernels include 
smoothing filters and edge detectors. OCT technology is relatively young, however, and 
extensive use has not yet been made of standard image processing techniques. 
6.2.2.1.2 Exponential correction 

A value set operation which is not linear and can not be implemented using a convolution : 
kernel is exponential correction. In the single-scattmng approximation, the detected OCT 
photodetector power fi:om a scattering medium attenuates vsdth depth according to ([28]; see also 
Chapter on Optical Coherence Microscopy): 

</^)^ oc Fiz) . 2Mtz] . (6.8) 
Here F(z) is a function of the focusing optics in the sample arm, m is the total attenuation 
coefficient of the sample (given by the sum of the absorption and scattering coefficients), and z is 
the depth into the sample. If the depth of focus of the sample arm optics is larger than several 
attenuation mean-free-paths (given by l/m) in the sample, then the function F(z) is relatively 
smooth over the available imaging depth and the attenuation may be considered to be dominated 
by the exponential temi. If this condition is not met (i.e. for imaging with high numerical 
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aperture), tlien the complete form of equation 6.8 must be taken into account. Equation 6.8 has 
been experimentally verified in model scattering media [28]. Thus, the reflectivity profile 
measured by OCT in a typical imaging situation is scaled by an exponential decay with depth. 
Because this decay is intuitively understood and expected, it is typically not corrected. It is 
possible, however, to correct the data such that a direct map of sample reflectivity is displayed. 
The analogous decay in ultrasoimd imaging is commonly corrected by varying the amplifier gain 
as a function of time by an amount corresponding to the decay ("time-gain compensation," or 
"TGC"). In principle, this approach could also be used in OCT by varying the gam with an 
exponential rise corresponding to the inverse of the expected exponential decay: e^^''^ where v is 
the speed of the depth scan. This approach can also be implemented after sampling by simply 
multiplying each A-scan point by point with an exponential rise, e^. This correction assumes, 
however, that the sample surface exactly corresponds to the first pixel of the Arscan, When not . 
true, this assimoiption will produce error, especially when the location of the tissue surface varies 
fi"om one A-scan to another. This error can be mitigated by first locating the sample surface, then 
applying the correction firom that location on: e^*^''^**^^, where z(0) is the location of the sample 
surface. Alternatively, it can be noted that the error amounts to a scaling error and the index of 
the surface location can be used to correct the scale. It should be noted that if the data has been 
logarithmically compressed, then the correction is simply a linear rise. It is clear that 
information is not added to the image by application of this type of correction, noise is sc2ded 
together with signal, and the deeper sample regions become extremely noisy. Therefore, it is a 
subjective matter whether exponential correction improves or worsens the image viewability, 
6.3 Signal Processing Approaches for Image Enhancement 
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6.3.1 Systems theory model for OCT 

In order to make use of signal and image processing concepts which are well known and 
ha\'c been extensively used in other medical imaging modalities, it is useful to describe OCT 
from a systems theory point of view. Fortunately, several authors have noted that since the low- 
coherence interferometer at the heart of OCT is basically an optical cross-correlator, it is 
straightforward to construct a simple transfer function model for OCT which treats the 
interaction of the low-coherence interferometer with the specimen as a linear time-invariant 
system [24, 43-45]. 

As the basis for the following sections on coherent signal processing in OCT, the model 
developed by Kulkami [24, 45] is here summarized. In this model, an optical wave with an 
electric field expressed using complex analytic signal representation as 
1e. - 2c, (c/ - z)e3q>[727iJfco ict - z)] is assumed to be incident on a Michelson interferometer, as 
illustrated in figure 6.15. We adopt a notation wherein variables ovetbarred by a tilde symbol (-) 
represent modulated quantities, whereas identical variables not so marked represent their 
complex envelopes. Thus, e/(cr-z) is the complex envelope of the electric field, is the 
central wave number of the source spectrum, and c is the free space speed of light The quantities 
/ and z represent the time and distance traveled by the wave, respectively. It is assxraied for the 
purpose of this model that the dispersion mismatch between the interferometer arms is 
negligible. By setting z = 0 at the beamsplitter, the field returning to the beam splitter from the 
reference mirror is given by 

V22;.(^2/^ ) = ^/(c/ - )exp[y 2;*o i^t - 2/^)] , (6-9) 
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where is the the reference arm optical path length. In the sample arm, the light backscattered 
from the sample with optical path length reaching the beam splitter is given by 

V2e; (^2/^ ) = V2e^ (cr - 21, )exp[y 2;*o ^ct - 21, )] , (6. 1 0) 

Here e/(c/-2/^)is the complex envelope of the backscattered wave. The fields returning from 
the reference and sample arms interfere at the detector to provide the superposition field 
= 2^ + The alternating component of the detector current is given by 

r(A/) oc(g; (c/^-2/Je;(cr-2/j) + (g:;(c/-2/,)e; (cr ~2/j}, (6.11) 

where the angular brackets indicate temporal integration over the response time of the detector 
and Al^2(lr'ls) is the round-trip path length difference between the interferometer arms. 

The source autocorrelation can be measured by monitoring the interferometdc signal when 
a perfect mirror is used as a specimen in the sample ami. We define the source interf^ometric 
autocorrelation function as: 

R'aiM)^(^i (cr-2/^)3}*(cf-2/j) = Rii(Al)oxpU27dco^O f (6.12) 

where Ra(Al)is the autocorrelation of the complex envelopes of the electric fields. The 
autocorrelation fijnction i?//(A/) is experimentally measured by demodulating the detected 
interferogram at the reference arm Doppler shift firequency, and recording the in-phase and 
quadrature components of this complex signal. According to the Wiener-BOiinchin theorem [46], 
the source spectrum is given by the Fourier transform of the interferometric autocorrelation 
function: 
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S,. ik) = Ji?;, (A/)exp[-y 2nM/]t/( A/) . (6.13) 
— « 

For an arbitrary specimen in the sample arm, the interferometric cross-correlation 
function is defined as 

^,(A/) ^ {e, ict - 2l^)e:ict - 2/,)) = i?,,(A/)exp(/2^oA/). , (6.14) 

where RisiAl) is the cross-correlation of the con:^)lex envelopes of the electric fields [46]. The 
cross-power spectrum S^^ik) also follows by analogy: 

CO 

S^(k)= J^^(A0expL/2it^(A0MA/). (6.15) 

6.3.1.1 Sample impulse response 

It is useful to represent the interaction between the specimen imder study and the incident 
sample ann light as a linear shift-invariant (LSI) system, characterized by a firequency dependent, 
transfer function [46-48]. The key insight this provides is that the transfer function H(lc) and its 
inverse Fourier transform, the sample impulse response h(z), are physically meaningful quantities 
describing the backscatter spectrum and electric field reflection coefficient profile of the sample, 
respectively. We denote the impulse response experienced by the electric field 2) by hiz) and 
that experienced by the complex envelope field by h(z) . This LSI model is valid when the 
contribution to fi"om multiply scattered light is not significant. 

The impulse response A(z) describes the actual locations and reflection coefficients of 
scattering sites within the sample, and convolves with the source electric field envelope to create 
the scattered electric field envelope: 
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eX-z) = 6, (-Z) (S>h\z). (6 J 6) 

Here ® represents the convolution operation, and the negative sign implies scattering in the 
negative direction (back-scattering). Taking the Fourier transform of both sides of equation 
(6.16) provides 

E,ik)^E^ik)H\k) (6.17) 

where £,(A:), E^^k), and if (X:) are the Fourier transforms of e^(r), e,. (z)and h(z), 
respectively. The assumption of shift invaiiance ensures that 

(ct - 2/ J = e, (ct - 2/, ) ® A * (-(c^ - 2/, )) . (6,18) 

Substituting equation (6.18) in the definition of the correlation functions (equations 6.12,6.14) 
provides 

RisW^RiiiAl)®h(Al). (6.19) 

The source autocorrelation function, Rfj(Al)and the cross-correlation function of the 
sample and reference arm electric fields, i?^ (A/) thus constitute the input and measured output, 
respectively, of an LSI system having the impulse response h(z). Therefore, the impulse response 
which describes the electric field-specimen interaction as a function of z is exactly the same as 
that which connects the auto-and cross-correlation functions of the interferometer as a function of 
the path-length difference -d/. Thus, this model provides access to understanding the fimdamental 
properties of the interaction of the sample wdth the sample arm electric fields by using simple 
correlation function measurements. 
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Equation 6,19 also leads directly to a simple, albeit naive approach for OCT image 
resolution improvement by deconvolution. Taking the Fourier transform of equation 6,19 and 
solving for the impulse response gives: 

where F'^ denotes the inverse Fourier transform, and S^(k) and 5// (A:) are the complex 
envelopes of the cross- and auto-power spectra, respectively. 

6.3.1.2 Depth-resolved spectral estimation 

The sample transfer function | H(k) \ vviiich may be obtained from the Fourier trasnsform 
of equation 6.19 describes the backscatter spectral characteristic of the sample, i.e. the ratio of 
the backscattered power spectrum to the spectrum which vvas incident However, since the 
coherence length in an OCT system is short, it is possible to use an analog of time-frequency 
analysis methods [49] to extract the backscatter characteristic with depth discrimination. This can 
be accomplished by limiting the detected transfer function data to the region of interest in the 

sample by digitally Mdndowing the auto- and cross-correlation data used to calculate \Hik)\ . 
This is the short-time Fourier transform [50] technique; wavelet transform [51] approaches may 
also be used. 

All information regarding the spatial and frequency dependence of scattering vwthin the 
sample is contadned in the complex space domain function h(z) and its Fourier transform Hik) . 

If the vnndowed region contains only a single scatterer, | H(k) |^ is straightforwardly interpreted 

as the wavenumber-dependent backscatter characteristic of that scatterer, denoted | C(^) | . In a 
medium containing many identical scatterers, the impulse response h(z) windowed to the region 
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of interest may be modeled as being a convolution of two new funclions: h(z) = b(z) ® c(z) , 
Here, b(z) describes the spatial distribution of scatterers along the sample axis z, and c(z) is the 
inverse Fourier transform of C(k). Under the assumption that b(z) is a white stochastic process 
(i.e., the scatterer locations are uncorrelated), the expectation value of | H(k) |^ is then given by: 

^{| P} = E{\ B(k) p} I C(k) M C(k) p , (6.21) 
where B(k) is the Fourier transform of b(z). Thus, the backscatter characteristic of the individual 
scatterers in a sample may be directly obtained within a user-selected region of the sample by 
appropriate Fourier-domain averaging of coherently detected windowed interferogram data. This 
analysis is readily extended to the case of a medium containing a heterogenous mixture of 
scatterers, each having its own backscatter characteristic. In this case, a similar signal processing 
algorithm produces an estimated spectrum corresponding to a weighted average of the individual 
backscatter spectra [45]. 

One final step allows for estimation of the actual backscattered spectrum of light returning 
firom the sample rather than the backscatter transfer characteristic of the scatterers. Since the 

spectrum of hght in the sample arm is defined as 5^^j(it) = ^j(A:)£*j(A:), 
usingE^ ik)^ Ef(k)H\k) (the modulated analog of equation 6.17) gives 

Sssif^) H 5^ Wr /Siiik), Similar to the previous paragraph, ensemble averaging must be used 
in order to average out apparent spectral features which are really due to the spatial structure of 
the scatterer distribution: 

Sss(k)^E{\SMh/S,,(k). (6,22) 
6.3.2 OCT image deconvolution 
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6.3.2.1 Straightforward approach 

The systems theory model described in the previous section predicts that the optical 

impulse response of tissue h{z) or h{z) is calculable if the complete cross-correlation sequence 
comprising the OCT signal is acquired with interferometric accuracy. The impulse response is 
interpreted as describing the actual locations and amplitudes of scattering sites wthin the sample 
arising from index of refraction inhomogeneities and particulate scatterers. 

An example of the application of Eq. (6.20) for direct deconvolution of undemodulated 
OCT A-scan data is provided in figxu-e.6,16. This data v^as acquired using a data acquisition 
system with interferometric calibration capable of capturing the crosscorrelation sequence with 
nanometer spatial resolution [24]. An interferogiam segment obtained with this systeni which 
includes several discrete reflections is plotted in the figure. Also plotted in the figure are the 
autocorrelatibn sequence from a mirror reflection, and the impulse response calculated using the 
modulated analog of Eq. (6.20). An increase in resolution by a fector Of >2 was obtained between 
the original interferogram and the calculated impulse response profile. TTie improvement 
obtained using this simple, no-cost algorithm is quite striking when executed on two-dimensional 
data sets, as illustrated in Fig. 6.17(a-b)- In this figure, digital deconvolution of magnitud&only 
demodulated A-scan data was used to improve image sharpness in the axizd (vertical) direction of 
a cross-sectional OCT image of a fresh onion specimen. 

The data used as the input for the deconvolution algorithm in figure 6.17a was acquned 
using an OCT system which, like most OCT systems constructed to date and most other 
biomedical imaging modalities (excluding ultrasound), records only the magnitude of the image 
data. A novel approach to signal deconvolution which tdkzs advantage of the complex nature of 
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OCT signals is to perform coherent deconvolution by supplying both the magnitude and phase of 
the demodulated A-scan data as complex inputs to the deconvolution equation. The advantage of 
this approach is illustrated m figure 6.18, which illustrates the capability of coherent 
deconvolution to extract real sample features from what would have been regarded as 
meaningless speckle in amplitude-only data. 

6.3.2.2 Iterative restoration methods 
The practical implementation of Fourier-domain deconvolution algorithms such as 
described by Eq. (6.20) inevitably leads to a reduction in image dynamic range, as evidenced in 
figure 6.17(b). This is because even small errors in the estimation of the auto- and cross- 
correlation spectra introduce noise in the resulting impulse response image due to the division 
operation inherent to deconvolution. In order to alleviate this problem, one study [48] has applied 
advanced constrained iterative deconvolution algorithms [52] to the OCT deconvolution 
problem. In these algorithms, improved restoration of the impulse response from the cross- 
correlation data is achieved by using prior knowledge of flie properties of the desired impulse 
response. This usage of a priori knowledge provides a consistent extrapolation of spectral 
components of the transfer frmction beyond the source bandwidth. The impulse response is then 
estimated by the method of successive approximations, hence no division operation leading to 
uistabilities and noise is involved. Therefore, the iterative algorithms offer the potential to 
achieve improvement in OCT resolution with a minimal increase in noise. In a test applying this 
algorithm to the same onion image data set displayed in Fig. 6.17(a), resolution enhancement 
similar to that displayed with the straightforward deconvolution algorithm was achieved, but 
with much smaller dynamic range loss (-2dB; see figure 6. 17(c)). 
6.3.2.2 CLEAN algorithm 
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One report has described promising results in deconvolution of OCT images using the 
CLEAN algorithm [44]. This algorithm, which is highly nonlinear and can only be described 
algoritlimically, begins with an estimate of a system's impulse response. It then searches a 
dataset for a peak, subtracts the impulse response from the data set at the location of the peak, 
and places a delta function at the location of the peak in the deconvolved dataset. Using a 
modified version of this algorithm, Schmitt demonstrated somewhat improved resolution and 
clarity in OCT images of phantoms and skin. 
6.3.3 Spectroscopic OCT 

Although most OCT systems developed to date have used light exclusively to probe the 
intemal miciostructure of sample specimens, the possibility of extracting spectroscopic 
information from OCT signals is particularly exciting since it may provide access to additional 
informatioa about the composition and ftmctional state of samples. Relatively little work has 
been published to date on this topic, primarily because the number of spectral regions in vMct 
sources suitable for OCT are quite limited, and because the spectral bandwidth of OCT sources 
are qiute narrow compared to light sources used in conventional spectroscopy. 

Efforts to date on spectroscopic OCT have concentrated in two areas. First, a few reports 
have appeared concerning spectral ratio imaging of OCT images acquired at two or more spectral 
bands [53, 54]. This approach can be implemented most elegantly by using a wavelengfli 
division multiplexer (WDM) to combine light from two or more sources into the source fiber, 
and then electronically distinguishing the resulting signals by their different Doppler shifts 
resulting from the reference arm scan. Of course, this approach is limited to a choice of 
wavelengths which can simultaneously be mode guided in a single-mode fiber, Ratiometric OCT 
imaging using a pair of sources at 1.3 and 1.5 microns (vrfiich are s^arated by approximately 
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one decade in water absorption coefficient, but have similar scattering coefficients in tissues) has 
been used to probe tlie water content of samples in three dimensions [53]. Combinations of otlier 
wavelength pairs have also been attempted in search of contrast in biological tissues [54], 

The second implementation of spectroscopic OCT is that described in section 6.3,1,2 
above, in which modifications of the source spectrum caused by the sample may be measured 
directly from Fourier-domain processing of cross-correlation interfeixjmetric data. The most 
successful application of this idea to date has been in Doppler OCT, in \^ich spatially resolved 
shifts in the sample spectrum due to sample motion are estimated from localized spectral shifts in 
the cross-correlation data [34, 35]. Details of the signal processing techniques used to extract this 
data and some preliminary applications are described in the Doppler OCT chapter. 

Even if sample motion is absent, however, the techniques described in section 6.3.1.2 may 
still be used to extract useful spectral information from samples [25, 37]. A simple example of 
depth-resolved spectroscopy using OCT is provided m figure 6.19. Here, the backscatter spectral 
characteristic | H (k) |^ of sample arm light reflected from the front and back surfaces (the latter 
having double-passed the filter) of a commercial interference filter are plotted. These quantities 
were obtained by applying the Fourier transform of eqmtion 6. 1 9 to an OCT A-scan of the filter, 
vrindowed to the vicinity of the glass-air reflection from the front and rear sides of the filter. It is 
interesting to note that this data illustrate quantitative depth-resolved reflectance spectroscopy 
demonstrating the equivalent of femtosecond temporal resolution in a simple, inexpensive 
system suitable for biomedical applications. 
6.4 Safety in the Clinical Environment 
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All important consideration in the implementation of OCT systems for medical diagnostic 
applications is tiie issue of operator and patient safety. All hospitals in the United States and most 
research funding agencies and foundations require Institutional Review Board (IRB) approval of 
all studies involving human tissues, including excised tissue samples. IRB approval procedures 
typically include the review of procedures for informed consent of patients who will be used as 
experimental subjects. Potential hazards to operators of OCT equipment and to patients fall into 
tliree categories: biohazards, medical device electrical safety, and optical radiation hazards. 
Procedures for avoidance of contamination and infection, as well as for electrical safety are well 
established in hospital environments, although ihey may not be as familiar to researchers vrfth 
physical science backgrounds. Biohazard avoidance primarily means utilization of proper 
procedures for. handling potentially infected tissues, as well as proper disinfection of probes and 
other devices which come into contact with patients or tissue samples. Electrical device safety, 
guidelines typically regulate the maximum current which a patient or operator may draw by 
touching any exposed part of a medical device, and are usually followed by including appropriate 
electrical isolation and shielding into the design of clinical OCT systems (see, for example, [55]). 
6.4.1 Optical radiation iiazards in OCT 

A potential operator and (primarily) patient safety concern which is unique to optical 
biomedical diagnostics devices is the potential for exposure to optical mdiation hazards. 
Although cw sources used for OCT are typically very weak compared to lasers used in physical 
science laboratories and even in other medical applications, the tight focussing of OCT probe 
beams which is required for high spatial image resolution does produce intensities approaching 
established optical exposure limits. A number of international bodies recommend human 
exposure limits for optical radiation; in the United States, one well-knovra set of guidelines for 
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optical radiation hazards are produced by the American National Standards Institute, ANSI 
Z136.I [56]. Unfortunately, these guidelines are specified for laser radiation exposure, and are 
also provided only for exposures to the eye and skin. Nonetheless, many analyses of OCT 
radiation safety have utilized these standards. The applicable ANSI standards for cw laser 
exposure to the eye and skin both recommend a maximum permissible exposure (MPE) 
expressed as a radiant exposure, which is a function of the exposure duration, and tabulated 
spectral correctiori factors. 



The following references are hereby incorporated herein in their 
entirety. 
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Multistep image dewarping of OCT images using 
Fermat's principle and mapping arrays 

All transformations will be carried out as back-transformations, meaning we determine 
for a given position (x,, y,) in the target image the position (x,, y,) in the raw data acquired 
by the frame grabber. In general we will use x for coordinates perpendicular to the A- 
scans and y for coordinates in line with the A-scans. For simplification, all coordinate 
systems have the origin in the center of the image. Pxxx is an abbreviation of (x^^x, yxxx) 
Dimensions of images are given in n pixels wide and m pixels high. 



Very often in imaging data is acquired as tomvs and colunms in a rectangular image. In 
many cases, this raw image has not the true aspect ratio, or more general, the resulting 
target image can have any shape. There are two ways to map the data from the raw image 
onto the target image. The forward transformations 



assign the point (xt, yO in the target image to where the pixel data at (x,, yO will 
correspond to. In general, (xt, yO won't be at an integer position, therefore the value at 
(xr, yr) has to be distributed on the neighboring pixels. TTie weights are not so ea^ to 
dertermine, especially if many raw pixels will be mapped to a small area in the target 
in^ge. The image has to be normalized against the pixel density. 
With the backward transformation 



the algorithm calculates for each pixel (xt,yt) in the target image the corresponding 
position (xr, y,) in the raw image. If diis position is not at a exact pfacel, fliere arc several 
ways to assign a value. The fastest ways would be the 'next neighbor', assigning the 
target pixel the value of the closest neighbor pixel of (xr. yr) in the raw image. Higher- 
precision can be obtained through bilinear iritMpolation between the four neighboring 
pixels. Other methods are trilinear or spline interpolation. 



yt=fytO(r.yr) 



(1.) 



(2.) 
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To do these transformation each time an image is acquired is, depending on the 
transformation, computational expensive. We are showing here a method that allows real- 
time image transformation using a mapping array and the next neighbor interpolation. 
The mapping array is an array of pointers, with the same number of rows and columns as 
the target image. If fxr and fyr are constant or seldom, the values of this array can be 
precalculated. The pointer at the position (xt, yt) will be assigned the address of the 
corresponding rounded pixel at the rounded position (Xr, yr). Once this has been done for 
all target pixels the image transformation can be done very quickly. To get the value for 
each target pixel the algorithm uses the corresponding pointer to access the pixel in the 
raw image (cf. Figure"?). Even complicated fxr and fyr do not slow down the imaging rate. 



B Polar image 

This technique is used in combination with the endoscopic probe. With this probe A- 
scans are taken in a radial feshion, with the probe constantly rotating. Therefore Xr and yr 
are rather polar coordinates: 

^ (3.a,b) 



with the radius fprobe of the probe and the imaging depth d. R and e are dimensionless. 
They can also be expressed in target coordinates 



(4.a,b) 



By combimng Eq. 3. and Eq. 4. the backward transformations can be obtained 



yr = 



^|xf+yf 



'probe 



'probe 



Si 
2 



Ulr. 

2n 



nt/2 



arcfan 



n. 



(5.a,b) 
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The image acquired by the frame grabber is called the raw image, with r as an index to 
define coordinates. Due to the sinusoidal motion of the reference arm mirror this image is 
deformed along the direction of the A-scan. Therefore the &st transformation necessary 
is from raw image coordinates 

The raw image is captured with nr pixels per A-scan and nw A-scans. In principle there 
would be n,pp pixel (peak to peak) available in an A-scan, therefore the duty cycle r\ is 
defoiedas 



11 = - 



(6.) 



The forward transformation from raw (Xr, yr) into intermediate (Xi, yi) coordinates is 
defined as 



y.(x.)=^sin 



2n 



2n 



(7.a,b) 



with the full peak to peak sunplitude mipp of the A-scan in the intennediate coordinates. 
Because tj < 100 % only part of the foil sinusoidal motion is visible, the intennediate 
image span m; pixel in depth, mipp can be calculated fiom 



m, 



(8.) 
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1 licrefore the back-transformation will be as follows 



n 



>in 



— arcsin 
7c y^ippj 



n, 



(9.a,b) 



0. ) Application on radially scanning probe 

For the combination of the polar image with the sinusoidal scanning Eq. 7. has to be 
inserted into Eq. 3*: 



mipp sin 



.fx 1 

Jin 

2 m; 



' probe ^ 2 



p/obe ^ 2 



2 sin I 
G(k,)=27c^ 



This we cian set equal with Eq. 4. and obtain 



X, = ^^arcsin 

71 




(10.a,b) 



2sin|-Ti 



+ d r. 



probe ^ 2 

d 



y, = ^arctara ^ 



J (ll.a,b) 
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When the handheld probe takes a B-scan, the scans emerges diverging from the final lens. 
The center of the image is aligned to be a focal length f away from this lens. The image 
scans a width w in the vertical center of the image and the scan depth d is measxared in the 
horizxDntal center. The A-scans are emerging radially from the focus, the pixel being 
narrower on top of the image than on the bottom. X| and yi are in kind of polar 
coordinates. In these coordinates, the angle <p and the distance L from the lens are defined 
in terms of the point Pi = (xu yi) 

<p(^;) = <Pmax;^ w!th<f^^arcsm(^] 

(12.a,b) 

nil o 

L is made dimensionless by dividing through mj. These <p and L can also be calculated for 
the target image: 



<p(>ff .y/)= arctanf-^^l with f, = 



(13.a,b) 



with L being dunensionless as well. <p and L can be set equal with Eq, 12. and these 
equations can be solved for (xj, y{): 
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(^f . ) = arctanf — ^ 



(14.a,b) 



This transformation can be combined with the correction of the nonlinear scanning of the 
reference arm (Eq. 9.): 



^r(^i./f)=— arcsin 



Tt 



(15.a,b) 



These equations are not dependant on the image content; therefore this transformation 
cem be realized as a mapping area, allowing for real-time imaging. 



If the mediimi that has been imaged is non-homogeneous in the index of refiraction, 
refraction occurs. Iii the following, we assnine that there are two smooth boimdaries 
between different media in the area imaged, with an index of refraction of no = 1 (air) on 
top, followed by two layers with the indices ni and ni. This can easily be expanded to 
several layers, in the expense of more computation time. In the current state the user 
defined the boundaries by a few points, which are than connected through a spline. 
The forward transformation would use Snell's law to calculate the target pixel given the 
raw data pixel. But for the back-transformation Fermat's principle has to be applied. It 
states that the light would always take the shortest patfi between die source and the target. 
In our case. 



XWID: <WO_0301 1764A2J_> 



wo 03/0 1 1 764 PCT/US02/24721 

68 



^^pXp ^'bi.y^v >^62.yz,2 ) = ^(^61. y^i. -^^,2./^^ 

has to be minimal. Li, L2, L3 are defined as 



(17.a,b,c) 



ybi and yb2 are functions of Xbi and Xb2, given by the user defined splines. Unfortunately, 
Xbi and Xb2 are unknown and had to be found through an optimization process to 
minimize L, which is computational intensive. Assuming that Xbi and Xb2 are not varying 
a lot between subsequent lines in the target image, this optimization can be simplified by 
taking the previous value as a seed and to look for the shortest path length if Xbi and Xb2 
are varied in steps of 0,1 pixel in the neighborhood of 0.5 pixel. When the first crossing 
point Pbi = (xbi, yb2) is found, <p can be computed: 

<p(Xi,i,yi,i)-arctanj^^^^j (18.) 

L and 9 from Eq. 16. and 18. can than be plugged into Eq. 12. to get the backward 
transformation: 

2<Pmax Vt-yb^J (19.a,b) 

y/(^r«yp^/,i.y/,i.^/,2.y/,2)=^/ -(^ +z^)m; 

and finally the complete back-transformation using Eq. 9, 



yrixtvy..)=^^ arctanf-^1 



(20.a,b) 



These transformations will be demonstrated in Figure//, using an image of the t^poral 
anterior camber angle as an example. 
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Correction of image distortions in OCT based on 

Fer mat's principle 



optical coherence tomography (OCT) is a relatively new tecMplbgjf/ yaiich is 
capable of mi^-on-scale resolution imaging noninvasively in living biological 
tissues, OCT shows rapid progress in the resolution, acqxiisition ra.te and possible 
applications. But one of the main advantages of OCT compared to ultrasound, non- 
contact imaging, also results in a ma or image distortion: refi:action at the air-tissue 
inter&ce. Additionally, applied scanning configurations can lead to deformed 
images, Botfi errors prevent accurate distance and angle ^measurements on OCT 
images. We describe a methodology for quantitative image correction in OCT 
which includes procedures for collection of non-telecentric scan pattens, as well as 
a novel approach for refitiction correction in layered media based on Fennat's 
principle. 



Introduction 
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Optical coherence tomography is a relatively new technology, which is capable of 
micron-scale resolution imaging noninvasively in living biological tissues. So far, the 
research focused on obtaining images in different applications (e.g. in ophthalmology, 
dermatology and gastroenterology), on resolution improvements real-time 
imaging, and on functional OCT like color Doppler qcT^^'^'*^^""*''** ^^^^^ or polarization 
sensitive OCT^^^^'^^-^^*^-^"^^^' * ^'^ K Meanwhile relatively little attention has been 
paid to image processing for quantitative image correction. Although it is relatively 
straightforward to use OCT to obtain accurate optical thickness measurements along a 
given axial scan line in a homogenous medium, there are many potential applications of 
OCT in which accurate two-dimensional reconstruction of complex features in layered 
media are required. One such application is in ophthalmic anterior segment 
biometry^^^^^"^ in which not only linear distances but also curvilinear 

surfaces, interfaces^ and enclosed areas must be accurately obtained. In this field, 
ultrasound biomicroscopy (UBM) has proven to be a value tool for the diagnosis of 
appositional angle closure, which is a risk factor for progressive trabecular damage, 
elevated intraocular pressure and acute angle-closure glaucoma^^^*^**' Since 
OCT is non-contact, imagmg the angle with OCT^'^*'^^"^***'-^^^^ greatly improves 
patient comfort, and allows for fast screening. An additional advantage is the substantial 
resolution increase from 50 to 10-15 pm. Unfortunately, the non-contact mode leads to 
strong image distortions due to refraction at the epithelium and and to lesser extend at the 
endothelium of the cornea. 

High lateral resolution at a convenient working distance of more the 10 mm complicates 
telecentric scanning*. The alternatively used diverging scan-pattem results in a 



SOCCID: <WO_030117e4A2_L> 



wo 03/01 1 764 PCT/US02/24721 

73 



second source of significant image distortions. Distortions due to non-paraxial beams 
could be reduced in other wide-angle imaging modalities like endoscopy 

Heircrty ci a!. 2001 ; Smith cl al. 1992) -xr ;r«o«;««(GfO"cnschiId 1997; Liu ct al. 2000) 

With the methods described in this letter we are able to reduce the maximum geometric 
errors due to non-telecentric scanning 10-fold down to 86 iim in a 3.77x4.00 mm^ image 
and achieve images corrected for refraction at multiple boundaries within the image 
applying Fermat's principle. 

Methods 

For display, P' = (x',/ ) in the acquired data, the raw image, has to be transformed into 
P = (x.y) target image. In principle, this can be done in forward (P = f(P*)) or 

backward (P' ^ F(P)) direction. For forward mapping the target position for a given data 
point is calculated This has a key disadvantage: Since the target position will most likely 
be between target pixels, sophisticated algorithms have to be applied to distribute its 
value onto the neighboring pixels to prevent dark spots and ambiguous assigned pixels, . . 
vAiich leads to a high computational exprase. Backward mapping avoids this 
disadvantage by mapping each target pixel to a location in the acquired image, then using 
simple interpolations to obtain its value. If the backward transfoimation is fixed, it can be 
implemented with lookup table to achieve real-time imaging^^^" In the raw 

image, x' and y' denote the coordinates across and along A-scans (single depth scans). To 
obtmn flie brightest possible images witii OCT, the field of view with a width w and 
depth d is centered a focal length f away from the lens on the optical axis. 
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Different scanning regimes can be differentiated, distinguished by the distance s between 
the pivot of the scanning beam and the final imaging lens with the focal length f (Fig. 
1 A). The position s - of the image of the pivot calculates to s*= (f • s)/(s - f), given by the 
thin lens equation. The scanning beams are either (i) converging (s>f, s'>f+d/2), (ii) 
parallel (telecentric, s=f), or (iii) diverging (s<f, s*<f-d/2). The condition given for s' 
avoids ambiguous assignments between raw and target image. P can also be defined in 
polar coordinates (<p,L), with the scanning angle (p and the distance L to a plane optically 
equidistance (EP) from the scaiming pivot. We have arbitrarily chosen this plane to 
intersect the optical axis in the lens (Fig. 1). In the general case and for a homogeneous 
sample (denoted by g and h), 9 and L are given by 

<pg(x,y) = arctan f 1 

Li; (x. y ) = s'-Vx^+(s^f-yF 
with the exertion of the teleoaitric scan (denoted by t), where s' goes to infini ty 

<Ph(x.y)=arctanf-^| 

\^-yJ . (2a,b) 

LUx.y) = n„Vx^+(f-yf 
The scanning angle to teach the outsides of the field of view at the focal plane is given by 

A=<(j.oj. (3) 

If the sample to be imaged is non-homogeneous, refraction occurs. Additionally, 
depending on the index of refraction n, the geometrical pathlength dififers from the optical 
pathlength. In the following, we assume k layers of different media in the area imaged. 



300CID: <WO_03011764A2LL> 



wo 03/01 1 764 PCT/US02/2472 1 

75 



with the indices of refraction of ni to nk, and air on top (no = 1). The boundaries Bk(x) are 
given as smooth functions. 

The forward transformation would use Snell's law to calculate the target pixel given the 
raw data pixel. But for the back-transformation Fermat's principle has to be applied. It 
states that the light would always take the shortest path between the source and the target. 
The pathlength can be divided into several pieces between the points Pj, where the beam 

crosses the boundaries, and the pathlength L^" in. air. <p is only depending on the first 
crossing point Pi. 

<p«-'(p,....p,.p)=<»(p,) 

L°«(P,....P,.P) = Lf{P,)^ E^'IP-PmI + n,|P,Pi 

Pi = (Xi,B,{Xj)) are a priory/^tially unknown. Fermat's principle states that the path 
length L of the light reaching P will be minimal. Assuming no foc^ spots, there is a 
unique solution for the Pi. 

In the rectangular array of acquired data, the horizontal position x' is linear with the scan 
angle 9', vMLo the equidistance plane is always a focal length away from the vertical 
center of ihe image: 

L'(x'.y') = f-y' 

Since cp*"'' = <p' and L^ * = L', the complete backward transfonnations x'= F^*(P,...,Pk.P) 
and y'= F^*(P,,..,P„,P) can be obtained from Eqs. (4), (5), assuming that <p^ = <f>^ : 
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F.S;'(Pi....Pk.P) = f-L*'(P,...,P,.P) 

F^'(Pi....P..P) = -^K'(P,) (6) 

The high speed OCT interferometer employed in this letter was based on a previously 
published design<'*°'"~«'^'- '^^«>, with a Fourier-domain rapid-scanning reference ann 
(RSOD) including a resonant scanner oscillating at fscan = 2 kHz. In the sample arm, we 
used a versatile lateral scanning hand held probe with a divergent scan (Fig. 1 Aiii) and a 
focal depth of 11 .4 ntm. This scanner was chosen for the best trade off between smallest 
focal spot sizes and large working distance for high patiait comfort The central image 
size was 3.77 mm wide and 4 mm deep (in air). A broadband (X,=1.32 fim, AX=68 nm, 
FWHM), high power (22 mW) semiconductor optical amplifier light source was used to 
illuminate flie interferometer. For improved SNR, we used an efficient, optical circulator- 
based approach with balanced, detection*""'"^*'"'- The images were preprocessed to 
remove the distortion form the nonlinear movement of the resonant scanner^^*^'^ " 

(36 fim maximum residual error). We imaged a test object and the temporal angle in 
the anterior segment of the eye of a healthy volunteer. Ail images were acquired at 
8 fi:ames/s to reduce motion artifacts. 

The transformations derived above were implemented using MatLab 5.2and applied 
offline to the acquired images in the following steps: First the geometric distortion was 
removed, therefore the first boundary could be distortion-firee defined semi-automaticaliy 
by user input of 4 to 6 points on the boundary and refined by active contours^^"""^ 
^^^K Second, after the correction of refraction at the first boundary, the second boundary 
was also distortion-fi-ee and could be defined. This scheme of defining boundaries and 
dewarping was continued until the image was completely dewarped. All intermediate and 
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the final target image always referred to the raw image data for minimum blurring due to 
the bilinear interpolation utilized. 

Since Fermat's principle has to be done for every pixel in the target image, this can be 
computational intensive if no previous knowledge is employed. We greatly improved the 
performance by assuming that the Pi for a given pixel in a row will only vary slightly 
from the PCs of the same pixel in the row before. Starting out with the previous Pi's, their 
Xi were first varied in steps of 0.5 pixels in a range of ±2.5. For the optimum path this was 
refined in a second iteration to a resolution of 0.1 pixels. 

Results 

As a first test object, we used a cover slip (975 jim, n2~1.52), placed on a piece of paper, 
with an Intralipid® drop (ni=L33) on top. Fig.l^Ai shows several distortions: (1) the 
boundaries of flie flat cover slip appeared bend, due to the geometric distortion of the 
diverging scarmer, and (2) under the drop the cover slip appeared to be bent down, both 
on the upper and lower surface, because the optical pathway to the bottom of the drop is 
longer than the physical. Maximum deviation from the flat sur&ce was 53 and.67 pm, but 
both effect partially compensated each other. (3) The cover slip showed up thicker than it 
physically was. Refraction was not obviously visible. After the correction, both cover slip 
surfaces were flat with a maximum error of 22 respectively 15 jun, the thickness of the 
cover slip was measured 963 ^un (Fig. Aii). Since the probe was hand-held, there is a 
remaining tilt due to non-normal positioning of the probe. Due to the highly positive 
curvature at the edges of the drop, two wedge-shaped areas are visible below, where the 
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• Snell's more sensitive lo non-smootliness of boundaries (if beam hits at locally step 
portion of b., the beam goes elsewhere), while Fermat's principle searches for absolute 
minimum 
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What is claimed is: 

1. A method for detemuning a condition of a tissue sample, comprising the steps of: 
generating at least one light beam; 

splitting the light beam into at least one measuring light beam and at least one reference 
light beam; 

directing the measuring light beam into a tissue sample; 

adjusting a relative optical path between the reference light beam and the measuring light 

beam; 

bringing the measuring light beam scattered back by the tissue sample into an 
interference relationship with the reference light beam; and 

processing the interferometric signal to determine a condition of the tissue sample. 

2. A systCTi for determining a condition of a tissue sample, comprising: 
light generating means for generating at least one U^t beam; 

beam splitter means for splitting.up said light beam into at least one measuring Ught 
beam and at least one reference light beam; 

adjusting means for adjusting a relative optical path between the measuring light beam 
and the reference light beam; 

directing means for directing the measuring light beam into the tissue sample; 

detection means for receiving Ught scattered back from the measuring light beam by the 
tissue sample; 

means for interferometrically superimposing the back-scattered measuring light beam and 
the reference light beam; and 

processing means for processing the interferometrically superimposed signal. 

3- Correction of multiple image distortions in real-time in OCT based on Fennat's 
principle and mapping arrays, 
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Figure 6.5. 
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Figure 6.9, 
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Table 6.1. Some useful spatial transformation (point-set operation) matrices. 
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